Genetic Drift in Genetic Algorithm Selection Schemes

A method for calculating genetic drift in terms of changing population fitness variance is presented. The method allows for an easy comparison of different selection schemes and exact analytical results are derived for traditional generational selection, steady-state selection with varying generation gap, a simple model of Eshelman's CHC algorithm, and evolution strategies. The effects of changing genetic drift on the convergence of a GA are demonstrated empirically

Annealing schedule from population dynamics

We introduce a dynamical annealing schedule for population-based optimization algorithms with mutation. On the basis of a statistical mechanics formulation of the population dynamics, the mutation rate adapts to a value maximizing expected rewards at each time step. Thereby, the mutation rate is eliminated as a free parameter from the algorithm.Comment: 6 pages RevTeX, 4 figures PostScript; to be published in Phys. Rev.

Genetic algorithm dynamics on a rugged landscape

The genetic algorithm is an optimization procedure motivated by biological evolution and is successfully applied to optimization problems in different areas. A statistical mechanics model for its dynamics is proposed based on the parent-child fitness correlation of the genetic operators, making it applicable to general fitness landscapes. It is compared to a recent model based on a maximum entropy ansatz. Finally it is applied to modeling the dynamics of a genetic algorithm on the rugged fitness landscape of the NK model.Comment: 10 pages RevTeX, 4 figures PostScrip

Error threshold in finite populations

A simple analytical framework to study the molecular quasispecies evolution of finite populations is proposed, in which the population is assumed to be a random combination of the constiyuent molecules in each generation,i.e., linkage disequilibrium at the population level is neglected. In particular, for the single-sharp-peak replication landscape we investigate the dependence of the error threshold on the population size and find that the replication accuracy at threshold increases linearly with the reciprocal of the population size for sufficiently large populations. Furthermore, in the deterministic limit our formulation yields the exact steady-state of the quasispecies model, indicating then the population composition is a random combination of the molecules.Comment: 14 pages and 4 figure

Evolutionary Games with Affine Fitness Functions: Applications to Cancer

We analyze the dynamics of evolutionary games in which fitness is defined as an affine function of the expected payoff and a constant contribution. The resulting inhomogeneous replicator equation has an homogeneous equivalent with modified payoffs. The affine terms also influence the stochastic dynamics of a two-strategy Moran model of a finite population. We then apply the affine fitness function in a model for tumor-normal cell interactions to determine which are the most successful tumor strategies. In order to analyze the dynamics of concurrent strategies within a tumor population, we extend the model to a three-strategy game involving distinct tumor cell types as well as normal cells. In this model, interaction with normal cells, in combination with an increased constant fitness, is the most effective way of establishing a population of tumor cells in normal tissue.Comment: The final publication is available at http://www.springerlink.com, http://dx.doi.org/10.1007/s13235-011-0029-

Authored Appendix to Cumulant Dynamics of a Population under Multiplicative Selection, Mutation, and Drift by M. Rattray and J. L. Shapiro

We revisit the classical population genetics model of a population evolving under multiplicative selection, mutation, and drift. The number of beneficial alleles in a multilocus system can be considered a trait under exponential selection. Equations of motion are derived for the cumulants of the trait distribution in the diffusion limit and under the assumption of linkage equilibrium. Because of the additive nature of cumulants, this reduces to the problem of determining equations of motion for the expected allele distribution cumulants at each locus. The cumulant equations form an infinite dimensional linear system and in an authored appendix Adam Prügel-Bennett provides a closed form expression for these equations. We derive approximate solutions which are shown to describe the dynamics well for a broad range of parameters. In particular, we introduce two approximate analytical solutions: (1) Perturbation theory is used to solve the dynamics for weak selection and arbitrary mutation rate. The resulting expansion for the system's eigenvalues reduces to the known diffusion theory results for the limiting cases with either mutation or selection absent. (2) For low mutation rates we observe a separation of time-scales between the slowest mode and the rest which allows us to develop an approximate analytical solution for the dominant slow mode. The solution is consistent with the perturbation theory result and provides a good approximation for much stronger selection intensities

The Dynamics of a Genetic Algorithm on a Model Hard Optimization Problem

A model of a hard optimization problem suggested in the literature is considered. The dynamics of a genetic algorithm (GA) using ranking selection, mutation and uniform crossover are completely modeled on this problem. These results are general and are valid for any symmetrical concave function of unitation. Full finite population effects are taken into account allowing a novel analytical comparison of roulette wheel and stochastic universal sampling. Closed form expressions are derived for the equilibrium population distribution of this model. The first passage time to move from a local to a global minimum in a two potential well landscape is calculated. A comparison is made with a stochastic hill climber and a GA without crossover. The GA with crossover is shown to perform orders of magnitude faster giving some insights into the nature of GA search and the crossover operator on these types of problem